Complex band structure with ultrasoft pseudopotentials: fee Ni and Ni nanowire 
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We generalize to magnetic transition metals the approach proposed by Choi and Ihm for cal- 
culating the complex band structure of periodic systems, a key ingredient for future calculations 
of conductivity of an open quantum system within the Landauer-Buttiker theory. The method is 
implemented with ultrasoft pseudopotentials and plane wave basis set in a DFT-LSDA ab-initio 
scheme. As a first example, we present the complex band structure of bulk fee Ni (which constitutes 
the tips of a Ni nanocontact) and monatomic Ni wire (the junction between two tips). Based on 
our results, we anticipate some features of the spin-dependent conductance in a Ni nanocontact. 
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INTRODUCTION 



Recent conductance data of nanocontacts and break 
junctions in magnetic transition metals such as Ni have 
shown interesting and partly unexpected results. While 
early conductance histograms for Ni at room tempera- 
ture in air appeared basically structureless 0, Oshima 
et al. 0, who worked in vacuum, at variable tempera- 
ture, and with the possibility of a magnetic field, found a 
minimal conductance step preferentially near 2 and 4 (in 
units of go — e^/Zi, the conductance quantum per spin) 
at RT and zero field, near 4 at 770 K and zero field, and 
near 3 (occasionally near 1) at RT with a field. Ono et 
al. reported again 2 for Ni in zero field, and 1 for 
Ni in a field. Break junction data by Yanson jj| show a 
minimal conductance step of about 3.2 in zero field. De- 
spite the poor consistency between these data, it is clear 
that the lowest conductance step for a Ni nanocontact is 
anomalously small in comparison with the large expected 
number oi s + d conducting channels, of the order of at 
least 8 for a monatomic contact. 

In a previous paper 0| we investigated the possibility 
that the smallness of the last conductance step observed 
for a Ni nanocontact may be caused by the fact that two 
tips are magnetized in different directions forming a mag- 
netization reversal (or domain wall) precisely in the con- 
tact region. By replacing the contact with a monatomic 
infinite tipless Ni wire we found that the number of con- 
ducting channels (energy bands crossing the Fermi level) 
can be significantly reduced from 8 (nonmagnetic wire) 
or 7 (ferromagnetic wire) to 2 as a magnetization rever- 
sal is built into the wire. The two remaining conduct- 
ing channels correspond to light s— electrons, while the 
heavy d— electrons are confined in the regions of positive 
and negative magnetization and cannot propagate along 
the wire due to the magnetization reversal. 

The real nanocontact consisting of two macroscopic 
tips and an atomic contact region is an open quantum 
system and has no artificial periodic conditions such as 
those in Ref. fs'l . The influence of the tips themselves (not 



considered in the previous model) will be in general im- 
portant for the nanocontact conductance. So, in order 
to make quantitative conductance predictions one will 
need to adopt a more detailed approach comprising two 
semi-infinite metals (either thick wires or semi-infinite 
crystals) connected by a neck (a short monatomic wire, 
or just a single atom). Several theoretical methods are 
available to study the transport properties of atomic scale 
conductors. Among these are methods based on nonequi- 
librium Green's functions combined with a localized basis 
set 0, 01 • The Landauer-Buttiker formula provides an- 
other simple way to calculate the conductance of quan- 
tum systems connected to two tips. The conductance is 
expressed in terms of the transmission coefficient at the 
Fermi level Ep- In order to compute the transmission 
coefficient one must solve a scattering problem for the 
open quantum system. The scattering approach has been 
applied by Lang and co-workers , and Tsukada and co- 
workersJolllOj using the jellium model for the electrodes. 
The layer KKR approach has been used to study the 
spin-dependent tunneling in magnetic tunnel junctions 
Lately Wortmann and co-workers have formu- 
lated the transfer matrix approach based on computation 
of single-electron Green function in the linearized aug- 
mented plane wave (LAPW) basis. We are interested in 
a formulation based on plane waves and ultrasoft pseu- 
dopotentials, a calculational technology which we are ap- 
plying to other problems involving transition metals |l3l | 
with great potential for the future. 

Recently Choi and Ihm presented a first-principles 
plane wave based solution to the scattering problem 
with real atomic contacts incorporating the Kleinman- 
Bylander-type pseudopotentials. We wish to generalize 
this scheme to deal with ultrasoft pseudopotentials 
which adequately describe the nuclei and core electrons of 
a transition metal such as Ni. The first step, on which we 
focus in this paper is the complex band structure of pe- 
riodic system, consisting of both propagating and expo- 
nentially decaying evanescent states. These generalized 
Bloch states will be needed to construct the scattering 
state deep inside the tips. The evanescent states play a 
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significant role in contact, surface or interface phenom- 
ena where there is a breaking of translational symmetry 
in one direction. 

We present here the complex band structure of bulk Ni 
(which gives information on the nanocontact tips) and of 
a monatomic Ni wire (the junction between the two tips). 
Although this is clearly only the first step towards a full 
theoretical description of the nanocontact conductivity, 
we shall see that some insight can be obtained already 
by examining the complex bands. The importance of 
complex band structure for conductance predictions has 
been recently emphasized for magnetic tunnel junctions 
[ifij i and molecular systems [T^ . 

METHOD 

We consider an open quantum system divided into 
three regions: left bulk contact {z < 0), scattering re- 
gion (0 < z < L) and right bulk contact (z > L) and 
we use a supercell geometry in the xy direction perpen- 
dicular to the wire. A wave propagating at energy E 
with k — (k_L , k) and incident from the left contact on 
the scattering region will form a scattering state 5* which 
due to the supercell geometry has a Bloch form in the xy 
direction: 
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^ik , R , 



^'k^(rj^,z;CT). (1) 



Different kj^ do not mix and can be considered separately. 
If there is no coupling between different spin polarizations 
in the Hamiltonian the same conclusion is valid also for 
the spin index a and we can omit both kx and a from 
the following formulas. 

Deep within the contacts the system becomes periodic 
also along the z direction and one can determine gener- 
alized Bloch states with complex k vector which obey to 
the periodic condition: 



-ipkir±,z + d) 
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■4^k{Y±,z), 



(2) 



where d is the length of the periodic unit cell deep within 
the contacts. In this region the scattering state ^ can be 
written as a combination of propagating and evanescent 
generalized Bloch states: 



in the left tip 
in the right tip 



k'eR 



where summation over k' G L {k' G R) includes the gen- 
eralized Bloch states in the left (right) contact at energy 
E which propagate or decay to the left (right). In the xy 
plane the wave functions ■0fc' obey the same Bloch condi- 
tion as in Eq. JQ) . Thus to compute the transmission and 
reflection coefficients {tkk' , ^kk' } one will need to find the 
complex band structure and generalized Bloch states tpk 
in the metallic contacts at the chosen energy E. 



We determine here the complex band structure within 
Density Functional Theory (DFT) describing the atoms 
with ultrasoft pseudopotentials |l5]. We consider the 
generalized Bloch states which satisfy the single-particle 
Kohn-Sham equation (atomic units — 2m = h = 1 
are used): 



ESliJk) 



Vnl 



(3) 



where Ves is the effective local potential (see Ref. 0, 
El), Vnl is the nonlocal part of the ultrasoft pseudopo- 
tential: 



Vnl - ^ D'^n\PL){Pn\ 



(4) 



Ivan 



constructed using the set of projector functions asso- 
ciated with atom / and S is the overlap operator: 
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\PL){f3i\ 



(5) 



The coefficients (7^„ are the integrals of the augmentation 
functions defined in Ref. 15|, il^]. It is convenient to 
rewrite Eq. (3) in the form: 
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From this equation it can be seen that the scattering 
problem in the ultrasoft pseudopotential case is not sub- 
stantially different from the norm conserving one (where 
E does not appear on the r.h.s.) since the energy E is 
fixed and is give n as an external parameter. 

As in Ref. [1J| we can write the solution of the integro- 
differential Eq. © as a linear combination: 
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(7) 



where ipn are linearly independent solutions of the homo- 
geneous equation: 



SV'„(r)= [-V2 + yeff(r)] ^„(r), 



(8) 



and 4'im is a particular solution of the inhomogeneous 
equation: 



S?A/m(r) = [-V2 + 14ff(r)] V^/™(r) 

+ ^e*'^^^^/3^(r-T^-Rj 



(9) 



Both t/jn and ^pIm are {x,y) periodic as in Eq. Q. Sum- 
mation over Im in Eq. (7) involves all the projectors in 
the unit cell and the coefficients cim,k are determined by: 

cim,k^J2^Di^n-EqU l[l3i{r-T')]*Mr)d'r. (10) 
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Note that the set of wave functions ipn is infinite. In prac- 
tice the expansion ■(/'n(r) = J2g^_ V'n(G±, 2:)e*(''^+*^^'''"^ 
over two-dimensional plane waves is used. If one consid- 
ers only N2D plane waves with G^^ < -Ecut the number of 
tpn becomes finite and equals to 2N2d- As in Ref.(l^] to 
find the functions ipn we discretize the unit cell along the 
z axis by dividing it into slices and replacing the Ves(j) 
in each slice by a z-independent potential. The ipn in the 
slices will be a linear combination of two exponentials 
with coefficients obtained by a matching procedure. 

The allowed values of k for a given energy E 
can be determined by imposing periodicity along z 
(Eq. (|2Jl) to the generalized Bloch function ipki^) — 
Y.G ^k{G±, z)e'('<^-i-+*^-i-)-''-L and to its z derivative: 



(11) 



(12) 



Inserting Eq. (7) into Eqs. (10)-(12) one can show that 
the last three equations are equivalent to the generalized 
eigenvalue problem: 



AX = e'^'^BX. 



(13) 



where A and B are matrices. We solve this prob- 
lem to obtain in general a complex k and coefficients 
X = |an,fc,c/m^fc| for the generalized Bloch state ipk at 
a given energy E and kj^. 

Our complex band structure calculations proceed in 
two steps. First we compute with a standard electronic 
structure code (PWSCF) the ground state electronic 
structure of the system (infinite bulk; and, separately, in- 
finite perfect nanowire) and we obtain the effective local 
potential 14ff(r) and the screened coefiicients I?^^ [2^ . 
In a second step, we use the potential T4ff(r) and the 
screened coefficients -D^„ to compute the complex band 
structure applying the scheme described above |2^ . 

RESULTS AND DISCUSSION 

As a first example, we studied the complex band struc- 
ture of ferromagnetic bulk fee Ni (Fig.l). The orienta- 
tion of the Ni tips in the actual break junction experi- 
ments is not under control, and we have chosen to plot 
the complex bands appropriate for a (001) orientation of 
the electrodes. We considered a tetragonal unit cell with 
two nickel atoms and the z axis parallel to the tetragonal 
axis. Our calculated optimal equilibrium lattice constant 
of ferromagnetic fee Ni is a = 3.42 A (expt. a — 3.52 A) 
and this value has been used for the band structure cal- 
culation. At this volume the magnetic moment per atom 
is 0.59 (expt. 0.61 ^b) H^]. In the scattering problem 
we took kx = for the purpose of examplification. Due 
to the symmetry of the Hamiltonian |25|, if at a fixed 




\mk 1.6 



FIG. 1: Complex band structure of fee Ni along [001] diree- 
tion at kx = for the majority (spin up) and minority (spin 
down) states. Tetragonal unit cell containing two atoms was 
used. Real bands (solid curves), complex bands with Re A; = 
(dotted curves) are plotted in the middle and left panels, re- 
spectively. The real projection of complex bands (open cir- 
cles) is plotted in the middle panel, and the corresponding 
imaginary projection (open circles) is plotted in the left or 
right panels. The band structures alog F-X direction in the 
fee Ni are shown on the insets. 



energy E a real A: is a solution of the problem, also —A: is 
a solution (corresponding to right- and left- propagating 
waves). If a complex fc is a solution also —k, k*, —k* are 
solutions (2^ (states with Im fc > (Im fc < 0) corre- 
spond to evanescent modes decaying along the positive 
(negative) z direction). We can therefore show the com- 
plex band structure just in the region Re fc > and 
Im fc > 0. 

Ni is a ferromagnetic metal and we plot in Fig. la and 
lb the bands corresponding to the majority spin (spin 
up) and minority spin (spin down) respectively. Figs, la 
and lb are divided into three panels. The propagating 
bands with Re fc ^ and Im fc = are plotted with con- 
tinuous lines in the central panel. Dotted lines in the left 
panel represent bands with purely imaginary fc (Re fc = 
and Im fc 7^ 0). The bands which correspond to complex 
fc with Re fc 7^ and Im fc 7^ are represented with 
open circles. They are shown with two branches which 
correspond to their projection in the real and imaginary 
planes. We chose arbitrarily to plot the projection on 
the imaginary plane in the left or in the right part of the 
figure in order to improve clarity. The band structure of 
fee Nickel corresponding to real fc obtained with PWSCF 
are shown in the inset of Fig. 1 along the F — X direc- 
tion. Since our cell has two atoms, the bands with real 
fc in Fig.l should coincide with the one in the inset after 
folding about the middle vertical line. We indeed found 
excellent agreement between the real band structures in 
the two calculations, with differences of the order of 0.005 
eV (invisible in the figure). 

In the band structure of fee Ni one can identify five 
narrow Zd bands and one broader 4s band. Along the 
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FIG. 2: Complex band structure of monatomic Ni wire for 
the majority (spin up) and minority (spin down) states. Real 
bands (solid curves), complex bands with Re k — 0, and vr/a 
(dotted curves) are plotted in the middle, left, and right pan- 
els, respectively. Each band is labeled by its main atomic 
character. 



(001) direction the s band mixes with the (1^2 band. The 
spin up d bands are fully occupied and only the s — dz^ 
band crosses the Fermi level. The down spin d bands 
are upshifted in energy by a calculated exchange split- 
ting of about 0.5 eV (exp. 0.33 eV). Therefore they are 
partially unoccupied and three d states in addition to 
the s — dz2 state crosses the Fermi level. The complex 
bands shown in Fig. 1 can be interpreted in terms of gen- 
eral theorems which describe their topological features 
|25|. It is known that whenever a band with real k has 
a local extremum there is a complex band with opposite 
curvature and the same symmetry which departs at the 
same Re k at right angles from the real plane. Complex 
k bands are parabolic or form loops joining real bands 
of the same symmetry. One such loop in the complex 
plane joins s — d^^ bands connecting the maximum at 
Re /c = 0.18 (27r/a) (very weakly pronounced) with the 
minimum at A: = 0.5 (27r/a). Comparing the complex 
band structures for different spin polarizations one can 
notice that there arc four propagating spin down bands 
at Ep while for spin up polarization the d bands cross the 
Fermi level only in the imaginary plane and are therefore 
evanescent (the smallest decay is Im k — 0.2 (27r/a)). 

In Fig. 2 we show the complex band structure for a 
monatomic Ni wire (a = 2.12 A). The supercell geometry 
and the other technical parameters used in this calcula- 
tion were already discussed in Ref. 0| . The bands at real 
k coincide (within 0.005 eV) with those presented there 
and calculated with the PWSCF code. Here all the oc- 
cupated real bands have maxima or minima at the zone 
border and the complex bands are therefore either purely 
imaginary (they are shown by dots in the left panel) or 
have Re fc = 7r/a (they are shown by dots in the right 
panel). Bands at imaginary k have a parabolic disper- 
sion except for the loop which joins dxz,yz bands. As 



in the case of bulk Ni at the Fermi energy there is only 
one real s — dz^ band of spin up polarization. All the d 
states are decaying and the smallest decay {dxz,yz states) 
is Imfc = 0.08 (27r/a). On the contrary 6 spin down bands 
with real k cross the Fermi level and the only decaying 
states crossing the Fermi energy have a much larger imag- 
inary part Im k = 0.35 (27r/a). 

The asymmetry between energy bands at Ep for differ- 
ent spin directions and the existence of evanescent states 
with long decay lengths (small decay parameters Imk) 
will play an important role in transport properties of the 
magnetic nanocontacts. One can see, for example, that 
if there is no spin reversal (two tips and the wire are 
magnetized in the same direction) all propagating states 
coming from the bulk Ni can get (though with some re- 
flection at the junctions) through the wire. On the other 
hand, if there is one spin reversal (situated either on one 
end of the wire or in the middle) the propagating d-states 
of one Ni tip become decaying in another Ni tip so that 
they will be completely blocked. The only open channels 
correspond to s-electrons (this conclusion was arrived at 
[3 considering a simple model of the Ni nanocontact) . 
And, finally, in the hypothetical case of two spin rever- 
sals (when the wire has the opposite direction of mag- 
netization with respect to both tips) the c?-electrons are 
only able to tunnel across the wire through the evanes- 
cent states. The latter have an exponential decaying form 
e~^/^ with the decay length £ = 7.96 a.u. {dxz,yz), 
5.31 a.u. (s — dz^), and 2.0 a.u. {d^y ^^-y^) which should 
be compared with the distance between atoms in the 
Ni wire a = 4.0 a.u. The shorter the wire, the larger 
the contribution of d-electrons to the total conductivity. 
This result is similar to the spin-dependent tunneling in 
magnetic tunnel junctions (see [HI Il6^) although in our 
case it is not so crucial since two propagating s-channels 
always exist. All these ingredients will become crucial 
when one will solve the scattering problem in order to 
find the transmission coefficient. It should be empha- 
sized that the number of propagating channels N give the 

2 

upper limit for conductivity, G = jj-N with unit trans- 
mission for all channels, while the scattering problem so- 
lution may produce very small transmission coefficients 
for some propagating channels. The channel number N 
of a thick wire will diverge proportional to the wire cross 
section, realizing the perfect conductance expected for a 
macroscopic case. 

CONCLUSIONS 

We are working to generalize to magnetic transition 
metals the approach proposed by Choi and Ihm for calcu- 
lating the ballistic conductance of open quantum systems 
within the Landauer-Buttiker theory. In this paper we 
focused on the complex band structure of periodic sys- 
tems, which is the key ingredient for the wave-function 
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matching between the scattering region and the metal- 
hc electrodes. The method has been implemented for 
the first time with ultrasoft pseudopotentials and plane 
waves in a DFT-LSDA ab-initio scheme. We showed that 
it can be applied to both the bulk solid and to the ide- 
alized monatomic nanowire. As a first application we 
calculated the complex band structure of fee Ni and of a 
ferromagnetic Ni monoatomic wire and discussed quali- 
tatively the properties of a nickel nanocontact which can 
be inferred from the complex band structure. Detailed 
calculations of conductance will be undertaken in forth- 
coming work. 
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